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Abstract 

Association schemes are combinatorial objects that allow us solve pro- 
blems in several branches of mathematics. They have been used in the study 
of permutation groups and graphs and also in the design of experiments, 
coding theory, partition designs etc. In this paper we show some techniques 
for computing properties of association schemes. The main framework arises 
from the fact that we can characterize completely the Bose-Mesner algebra 
in terms of a zero-dimensional ideal. A Grobner basis of this ideal can be 
easily derived without the use of Buchberger algorithm in an efficient way. 
From this statement, some nice relations arise between the treatment of zero- 
dimensional ideals by reordering techniques (FGLM techniques) and some 
properties of the schemes such as P-polynomiality, and minimal generators of 
the algebra. 



1 Overview and background. 

In this paper the author is interested in the development of the connection between 
the algebraic properties of the Bose-Mesner algebra associated to an association 
scheme and the combinatorial properties of the last one. This investigation is 
mainly motivated by the fact that computing the eigenvalues of an association 
scheme is in general difficult (even numerically). We are interested in applications 
of the theory of association schemes to the areas of coding theory and design theory 
[@, [|, @, |[ H^, therefore, we mainly treat commutative symmetric association 
schemes. Anyway we do not discard the non commutative case and some directions 
are showed in the conclusions. Main tools the author proposes to show the above 
connection is using some reordering techniques from zero-dimensional ideals called 
FGLM techniques. The main objective is pointing the connection between the 
computer algebra treatment of zero-dimensional ideals and association schemes. 
In fact much about metric properties (P-polynomial properties) and Q-properties 
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can be shown just by considering an appropriate set of equations. The paper is 
organized as follows: in the first section basic facts on association schemes as well 
as FGLM techniques are dealt with. We mainly follow the definition of an association 
scheme given in m. For an account on association schemes see for instance |L5l 



section we build a correspondence between properties found by FGLM techniques 
and combinatorial properties of association schemes. Most of the new material 
is placed in this last part, even though some approaches to known combinatorial 
results in section 2 are also new. Finally in the conclusions section we show some 
open questions and lines of research. The approach in this paper is a computational 
one, as a final remark to this introduction we point that some other computational 
approaches can be found in motivated by the works in [jl7[ |l9| . 

1.1 Association schemes. 

Throughout this paper we will follow the definition of association scheme given by 
E. Bannai and T. Ito (l): 

Definition 1. An association scheme with d classes is a pair S = (X,{Ri}f =0 ), 
given by a finite set X and a set of relations {Ri}f =a on X , satisfying the following 
rules: 

1. Rq = {(x,x) \ x E X} (the diagonal relation) 

2. {Ri}f =0 is a partition on X x X. 

3. Vi € {0, . . . d} 3j G {0, . . . d} such that R\ = R Jt where 

R\ = {(»,*) I (af.y) e Ri} 

4. For each election ofi,j, k S {0, . .. d}, the number: 



is constant for all (x, y) € Rk 

We can rewrite the above conditions as matrix relations. Consider the set of 
square matrices of order v — \X\ given by: 





\{zeX\ (x,z) e Ri (z,y)eRj}\ 



i = l,... 



,d Di = [Di(x,y)]o<x, v <v Di{x,y) = 



1 if (x, y) e Ri 



elsewhere 



Therefore, the conditions 1 — 4 above are equivalent to: 



1 '. Dq = Id (identity matrix) 



2'. X)fe=o Dk — Jt where J is the matrix where all entries are 1. 
3'. Vi S {0, . . . d} 3j e {0, . . . d} such that D\ = D 3 

4'. A-^=Eto4^ 
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We say that the scheme is commutative if p\^ = p k ^ Vi,j,k G {0, . . .d}. We 
say the scheme is symmetric if R\ = Ri. The set of matrices {Di}f =0 is the 
generating set of a semisimple algebra B over C called Bose-Mesner algebra or 
adjacency algebra. All matrices in the set are linearly independent by condition 
2', hence B has dimension d + 1. If the scheme is commutative, it is clear that 
DiDj = DjDi, and therefore the algebra is commutative. From now on, for the sake 
of brevity, we will write scheme for a commutative symmetric association scheme. 

Let / be the identity matrix and J the v x v matrix where all the entries are 1. 
It is well known (see [Q) that B is a diagonalizable algebrafjand there is a unique 
set of primitive idempotents {Eq = jy\J, E\i . . . , Ed}. 

Definition 2. The matrix representing the change of base from the base of the 
adjacency matrices to the idempotent matrices is called character table?] and it is 
denoted by: 

Note that the entries on the i th row are the eigenvalues (counting multiplicities) 
of the matrix D^. Also, since we are considering symmetric association schemes, 
the Di's are symmetric and the eigenvalues are real. If we consider the algebra U 
generated by {E , Ex, . . . , Ed} endowed with Hadamard product over the complex 
numbers, we have that the set of primitive idempotents is the set {Di}f =Q (see 
and: 

d 

Ej = fcCflA q = fornix = \x\p- 1 (2) 

2 = 

The qi(j)'s are called second eigenvalues. 

1.2 Grobner bases techniques. 

In this paper we use some well-known facts on Grobner bases theory, for background 
material on this theory see ^ or |24). Main topics we need is the relation between 
diagonalization and solving O-dimensional ideals and some clues on FGLM results. 
Treatment of finite generated algebras can be found in || |l4], 

1.2.1 Solving O-dimensional ideals and Grobner bases. 

First we show some relations of the problem of solving a O-dimensional system of 
equations with an eigenvalue problem. This material can be found in |2^, EjJ. Indeed 
this fact discloses the relation between the systems proposed in the next section and 
the problem of diagonalizing the Bose-Mesner algebra. Roughly speaking it works 
as follows: 

We are given an ideal X generated by T, where T is a reduced Grobner basis 
for T in some monomial ordering and V(J-) (the variety generated by J 7 ) is zero- 
dimensional. As usual lt(f) will denote the leading term of / and the set of terms 
will be denoted by: 

T = {4°...a# |* ,...,*«ieNo} (3) 

1 This arises from the fact that {D{\ is a set of commuting matrices. 
2 The justification for this name is given in |25|. 
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We define the normal set as: 

N = {t e T I $f e T such that lt(f)\t} (4) 

N is a basis of C[xo, . . . , xj] modi and the cardinality of N is just the number 
of roots of T . We can describe the effect of multiplying an arbitrary element by 
a fixed / G C[xo, ■ ■ ■ , x<j\ modulo I just by using the images of multiplying / by 
each term in N. Hence, if we denote by nf (/) the normal form with respect T , we 
conclude that: 

nf (/ • Xj) = Qjj(/)ar, i = 0, ...d (5) 

As usual the complex matrix A(f) = [a,ij(f)]f j =0 is the multiplication matrix of 
/. The multiplication tables for the elements in N are particularly important since 
from them we can rebuild any multiplication table and clearly A(xi) — \Pij]j ^—q- 
The following result is a slightly modified version of a theorem in |22| . 

Theorem 1 (Moller, Stetter). Let T a Grobner basis, Mq, . . . , Md the multipli- 
cation matrices for xq, . . . , Xd defined as above. Let U the unitary matrix such 
that U*MiU = Fi is diagonal for each i e {0, 1, ... , d}. If we denote the diag- 
onal entries of the matrices Fi as (uq, . . . , u d ) then the points on V(J-) counting 
multiplicity are: 

Zj = (x 3 07 ...,x 3 d ) = (u°, «},..., uj) (6) 
i.e. the roots of T are the eigenvalues of the Mi. 

Proof. For a proof of the result see [jz^] . □ 

Therefore, as usual, we will reduce the problem of diagonalizing the algebra B to 
solve a system of equations in d indeterminates over C. 

1.2.2 FGLM results. 

Finally in this section we point here briefly to two properties on Grobner bases for 
0-dimensional ideals. An introductory text to these ideas can be seen in Q or 
[ p4[ , the techniques and results exposed here can be found in [[13], We say 
that an ideal / 6 R[xo, . . . ,x n ] is in generic position with respect n if, for each 
(drj, ■ ■ ■ , o n ) and (b , . . . , b n ) distinct zeros of / in R", a„ ^ b n . We say that a 
change of coordinates is a generic change of coordinates for an ideal / if it is in 
generic position for one of the variables. 

Proposition 1 (Gianni, Mora). Let I be a zero-dim. ideal in M[xq, . . . ,x n ], there 
is a generic change of coordinates for I, and there are g , . . . ,g n e R[x„] such that 
if i < n then degree(gi) < degree(g n ) = d where d is the number of distinct zeros 
of I and g n is squarefree, and finally 

G = {x - g ,x 1 -gi,... ,x n -i - g n -i,g n } 

a reduced Grobner basis of I w.r.t. the lexicographical ordering such that xo > 

X\ > ■ ■ ■ > x n 

Proof. See the proof of the first assumption in proposition 1.6 in [Q. □ 
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Proposition 2. There is a probabilistic algorithmic method of computing such a 
generic change of coordinates. 

Proof. For a proof see Jl^], in fact, that paper is devoted to compute such an 
expresion. □ 

We will show the combinatorial significance of these facts in the case of the Bose- 
Mesner algebra later in this paper. 



2 Ideals related to association schemes. 

In this section we introduce some systems of equations having as solutions the first 
and second eigenvalues of the scheme. We will interpret some results both in terms 
of computer algebra and combinatorics, showing the great resemblance between 
both topics. 



2.1 First eigenvalues. Interpretation of the intersection alge- 
bra. 

It is a well known fact in combinatorics that the first eigenvalues of the scheme 
defined above are the nonzero solutions of the system given by: 

d 

XjXj - ^PjjXk = 0<i<j<d (7) 

fc=0 

This fact can be proved in many ways, we will sketch a proof in a Grobner basis 
flavor. For a complete and rigorous proof see here we will be more concerned 
with showing the parallel noted above. 

Let us consider the set of polynomials (as usual x = 1, i.e. the identity relation): 

d 

T = {xiXj - p° - ^2 Pij x k}i<i<j<d (8) 
fe=i 

From now on we will call equations in (7) structure equations and the polynomials 
in T structure polynomials. There are some nice properties on T the ideal generated 
by T : 

1. The variety V(X) is just the set of nonzero solutions of the equations in (7). 

2. T is a reduced Grobner basis for (T) = X in the total degree monomial 
ordering. 

3. T is a radical ideal, ie. the variety V(I) is reduced. 

We recall the above section and build up the multiplication matrices A(f) — 
[dij (f)]f j— for / G T. The multiplication tables for the elements in the corre- 
sponding normal set N = {1, x\, . . . , Xd} are the matrices given by M l , < i < d 
where: 

M i = {m\ m ) d l m=0 , where m\ m =pu ( 9 ) 
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These matrices are known in combinatorics as intersection matrices and it is well 
known that they generate an algebra isomorphic to the Bose-Mesner algebra called 
intersection algebra. Solving the equations in (7) is the same as diagonalising the 
intersection algebra and, since they are isomorphic, is the same as diagonalising 
the Bose Mesner algebra. This is the main fact that offers us a bridge between 
properties of these algebras and properties of the ideal X. 

2.2 On the other hand ... 

Another well known fact is that the solutions of the d quadratic equations: 

d d 

x k - J2 p oj Xk ~ Pij x i x i = 0<k<d (10) 

3=1 i,3=l 

is the variety of idempotents or projectors of the algebra. I.e. the solutions are 
the second eigenvalues of the scheme. The main trade off of this system is that 
the polynomials involved do not form a Grobner basis of the variety and in general, 
computing Grobner basis is a non-efficient problem, but once achieved a Grobner 
basis of it, all computations are carried out in the same way. Some tricks to deal 
with these equations are sumarized in the reference [Q. We will consider only the 
first eigenvalues from now on. 



3 Combinatorial properties achieved by reordering 
techniques. 

One of the main facts when dealing with association schemes is to know whether 
they are P-polynomial or not. We say that a scheme S = (X, {Ri}f =0 ) is P- 
polynomial if we can reorder the relations Ri so that the corresponding Di is a 
polynomial in Di with degree i and p(xx) is a polynomial in X\ whose roots are 
the eigenvalues of D\. 

Clearly, if S is a P-polynomial scheme we have that there is a monomial ordering 
such that the set: 

{^0 - 1,^(^1)1^2 ~P2(x il ), ...X id -Pd(xij} (11) 

is a reduced Grobner basis for T for a pure lexicographical order where Xi 1 <Xi , j = 
2, . . . , d. It follows that next proposition holds: 

Proposition 3. Let S = (X, {Ri}f =0 ) be an association scheme and X its associ- 
ated ideal. S is a P-polynomial scheme if and only if there is an pure lexicographical 
ordering for the variables Xi j with x io = xq, x^ < Xi j ,j^io such that 

{x Q - l,p{x il ) } x i2 -p 2 (x il ), ...x id -Pd(xij} (12) 

is a Grobner basis fori for that ordering, and the following conditions on the degree 
of the polynomials p,pi, . . . ,pd are satisfied: 

degree(pi) = i Vi = 2, . . . , d degree(p) = d (13) 
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Note that for an association scheme to be P-polynomial its Bose-Mesner algebra 
is the minimal algebra containing one of the relations but the converse is not true in 
general, see example 19 in |2l) . The reason for this is the conditions on the degrees 
of the polynomials. Note also that the resulting univariate polynomial is squarefree 
and it is just the characteristic polynomial of the relation that "spans" the whole 
algebra^. The main profit we get from proposition above is an algorithmic way of 
checking the P-polynomiality of an association scheme: 

Algorithm 1 (Checking P-polynomiality). 

• Input T Structure equations of S, 

• For i from 1 to d do 

Compute a reduced Grobner basis for (T) for a pure lexicographical monomial 

ordering where Xi < Xj j ^= i 

Check the conditions on proposition above. 

If the conditions are fulfilled then Stop and 

Return The scheme is P-polynomial and the Grobner basis computed fi 
od 

Stop Return The scheme is not P-polynomial 



Complexity of the algorithm is placed in the Grobner basis computation, but 
since we are dealing O-dimensional ideals and T is already a Grobner basis, these 
computations can be done by FGLM techniques which are effective and efficient. 
Clearly the complexity is at most the same as d FGLM Grobner basis computations. 
The algorithm can be easily modified to check whether the scheme is P-polynomial 
with respect to more than one relation or to check the minimal number of relations 
spaning the whole algebra. This last fact is important in coding theory for computing 
the Lloyd polynomials and metric properties of the scheme [Q. In that paper can 
be found many examples of how to apply this algorithm, we recall just two of them: 

Example 1. Let X = Z m ,m G N and a radix r (r < m) . The subgroup of7L* m 
generated by r and —1 acts on X by multiplication with orbits Ok- The following 
relation defines an association scheme (2-orbit association scheme) on X : 

xR k y Oi-i/eOt 

When m = 9,r = 2 the orbits are: O = {0}, 0\ = {-1,1 - 2,2,-4,4}£> 2 = 
{—3,3} T is given by: 

T = {x2 2 — 6 — 6x1 — 3x2, xq — 1, x\ 2 — x\ — 2, x\X2 — 2x2} 

The above algorithm return us the following Grobner basis: 

{—X2 2 + 6 + 6x1 + 3x2, xo — 1, — 18x2 — 3 X2 2 + X2 3 } 

Which allows us to compute the character table by solving the equation —I8X2 — 
3x2 2 + x 2 3 = and evaluating the solutions in the other expressions: 

/ 1 2 6 

CharTab = 1 2-3 

V 1 -1 

3 This means that any element in the algebra can be expressed by a formula of the generator 
and J, using the addition, the ordinary multiplication and the Hadamard product. 
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The following example is not P-polynomial: 

Example 2. Consider the scheme defined as above given by m — 8, r — 3, the 
orbits are: O = {0}, d = {±1,±3}, 2 = {±2,}, 3 = {4} The associated 
ideal X is generated by: 

{x x 2 - 4a; 2 - 4x 3 - 4, x 3 2 - 1, Xq - 1, x x x 2 - 2s b 
x 3 -xi, x 2 2 - 2x 3 - 2, x 2 x 3 - x 2 } 

If we compute the reduced Grobner basis with respect plex. orders we have: 

1. X\ < x 2 , x 3 

4^3 + 4^2 — x\ 2 + 4, 2x2 2 + 4x2 — %i 2 , x\X2 — 2xi, — 16xi + xi 3 , xo — 1 

2. x 2 < xi, x 3 

x - 1,X! 2 - 4x 2 - 2x 2 2 , xix 2 - 2xi, -x 2 2 + 2x 3 + 2, -4x 2 + x 2 3 

3. x 3 < xi, x 2 

xi 2 - 4x 2 - 4x 3 - 4, xix 2 — 2xi, xix 3 — xi, x 2 2 - 2x 3 - 2, x 2 x 3 — x 2 , 
x 3 2 - l,x - 1 

None of the three cases above fulfill the conditions on the proposition, and therefore 
is not P-polynomial. 

In fact, as we said above, a small modification on the algorithm allows us to 
compute the minimal number of relations spanning the algebra: 

Example 3. Following the setting in the example above, if we choose one ordering 
(for example x\ < x 2 , x 3 ) we have the following expression in two variables of the 
relations: 

x = l,xi = xi,x 2 = x 2 , x 3 = ^(xi 2 - 4x 2 + 4) 

Finally, propositions 1 and 2 allows us find an algorithmic procedure to find a 
non-negative integer matrix A in the Bose Mesner algebra such that any element in 
the algebra is expressible in terms of A, i.e. spans the whole algebra]]. Clearly, this 
happens if and only if there is a change of coordinates such that A is a resulting 
coordinate, and I is in generic position w.r.t. A. Therefore lemma 1.1.1 in [|l7| 
becomes obvious in a Grobner basis setting also as they recall in their paper from 
a linear algebra approach. This matrix is of great importance when we compute 
partition designs on the association scheme. A modification of algorithm 5.4 in [Q 
allows us find such an element A, the idea is making changes of coordinates until the 
ideal is in generic position, this is achieved by proposition 2. For the understanding 
of the modifications we made in the algorithm one must take in mind that our 
ideal is already radical, so we do not need to compute the radical, and also that 
we are looking for an non-negative integer change of coordinates which imposes no 
restriction on the running of the algorithm (For a proof see 

4 The meaning of this matrix A in terms of graph theory can be found in [p. 71] 
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Algorithm 2 (Finding the matrix A). 



• 


Input T Structure equations of S, 


• 


Compute a reduced Grobner basis Q for {T) for a pure lexicographical mono- 




mial ordering where Xd < x 3 . < j < d 


• 


Repeat 




Let {g} =gnR[x d ]. 




Let i the maximal index s.t. for j < i, Xj is expressible in terms ofxd- 




If i < d — 1 then 




Choose a random positive integer c 




Make the change of coordinates 5 




x i l— * x i < j < d 




CXi + X d H-> x d 




Compute a reduced Grobner basis Q with the same ordering in the new vari- 




ables, fi 




Until i = d — 1 Stop 


• 


Return The Grobner basis and the change of coordinates to achieve the last 




x d - 



Readers interested in the algorithm are referred to the paper [jl2| for a more 
detailed description of why it works. We just show how to deal with example 2 
above: 



Example 4. Taking the change of coordinates given by: 

x i | — > Hi < i < 3 

2^3 + %1 >-> V3 

The structure equations in the new base become: 

vl - 4y 2 - 4(y 3 - yx) ~ 4, (y 3 ~ z/i) 2 -l,y - 1, y x y 2 - 2yi 
yi(V3 - Vi) - 2/i, 2/1 - 2(2/3 - 2/1) - 2, 2/2(2/3 - 2/i) - 2/2 
/4nc/ a Grobner basis for a plex. ordering where j/3 < 2/1, 2/2 is given by: 

12yi - j/3 3 + 3y 3 2 + t/ 3 - 3, 27 + 24y 2 - 3y 2 + 25y 3 - y 3 3 , 
15-16yl + 2y 3 -2y 3 3 +yl,x a -l 

Hence the ideal is in generic position and the algebra is expressible in terms of 
the matrix given by the sum of the relation D\ and relation D3. 

4 Conclusions 

Properties of symmetric commutative association schemes can be derived from 
their associated algebras via ideal computations and Grobner basis of the ideals 
X associated with then allow us to describe a dictionary between combinatorial 
properties and algebraic treatment of zero-dimensional ideals (see appendix). We 
have proved that these ideals are radical, and therefore any polynomial condition 
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on the points V(I) can be seen as a membership problem which is answered by 
computing the normal form. All computations and algorithms proposed are effective 
and efficient since we are dealing with zero-dimensional Grobner basis using FGLM 
techniques |To| on the ideal. It is not surprise one can use Grobner methods here 
since we are dealing with finite dimensional algebras, and the linear algebra in the 
background is indeed that those researchers working on association schemes use. 
The main contribution of this paper is show this connection and also the algorithms 
on the paper give some extra information about the association scheme for free. 
Algoritm 1 does not only check P-polynomiality but gives us the minimal number 
of relations generating the algebra and the expresion of the Lloyd's polynomials and 
the matrix A resulting from algoritm 2 gives us a valuable tool for checking partition 
designs on the scheme. 

Moreover, future trends of work point towards extending this setting to the 
non-commutative case using the FGLM techniques in We also would like to 
study the relationship pointed by the Galois correspondence studied in and the 
corresponding ideals. Finally a new direction of work is develop a similar approach 
for dealing with proper and adjacency polynomials in the theory of graph spectra 

0- 

5 Appendix. 

This paper allow us compose a small dictionary translating combinatorial properties 
of the association scheme and properties and algorithms applied to is associated 
O-d imensional ideal. It can be set down as follows: 

Association Schemes Algebraic theory of finite algebras. 

S=(X,{R i }? =0 ) C[x ,...,x d ]/(F) 
Character table of <S Points on V{(T)) 

P-Polynomiality Algorithm 1. 

Minimal fl of relations expressing the scheme Minimal fl of generator for the algebra 

Algorithm 1 modified. 
Single matrix for expressing the scheme General change of coordinates for (JF) 

Algorithm 2. 

Metric properties 
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